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Abstract 

In this paper, we address the design of codes which achieve modulation diversity in block fading single- 
input single-output (SISO) channels with signal quantization at receiver and low-complexity decoding. With an 
unquantized receiver, coding based on algebraic rotations is known to achieve modulation coding diversity. On 

> 

' the other hand, with a quantized receiver, algebraic rotations may not guarantee diversity. Through analysis, 

in : 

we propose specific rotations which result in the codewords having equidistant component-wise projections. 
T^lj- ' We show that the proposed coding scheme achieves maximum modulation diversity with a low-complexity 

o 

minimum distance decoder and perfect channel knowledge. Relaxing the perfect channel knowledge assumption 
we propose a novel training/estimation and receiver control technique to estimate the channel. We show that our 
\ coding/training/estimation scheme and minimum distance decoding achieve an error probability performance similar 

J;? ' 

■ to that achieved with perfect channel knowledge. 
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I. Introduction 

In practical communication receivers, the analog received signal is quantized into a finite number of bits 
for further digital baseband processing. With increasing bandwidth requirements of modem communication 
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systems, analog-to-digital converters (ADC) are required to operate at high frequencies. However, at 
high operating frequencies, the precision of ADC's is limited [[B. Limited precision generally leads to 
high quantization noise, which degrades performance. In case of fading channels, floors in the bit error 
performance have been reported, and it seems difficult to avoid this behavior [|2l[|3l. On the other hand, 
channel capacity results show that even with 2-bit quantizers, the capacity of a quantized output channel 
is not far from that of a channel with unquantized output [@1[|5||. Therefore, there appears to be a gap 
between the theoretical limits of communication with quantized receivers, and the current state of art. 

In communication systems with fading, an important performance metric is the reliability of reception. 
For single antenna fading scenarios, modulation diversity is a well known signal space diversity technique 
to improve the reliability /diversity of reception [|6][E]. However, with a quantized receiver, this coding 
alone does not guarantee improvement in diversity. 

In this paper, we propose 2-dimensional constellations rotated by an angle 6 which can achieve full 
modulation diversity with a quantized receiver. With a quantized receiver, the maximum likelihood (ML) 
decoder is not the usual minimum distance decoder, and would be much more complex to implement. We 
therefore assume a minimum distance decoder operating on the quantized channel outputs. We observe 
that, with a quantized receiver, /) for a given rate of information transmission in bits per channel use, 
there is a minimum requirement on the number of quantization bits, without which floorjl] appear in the 
error probability performance, ii) there is only a small subset of admissible rotation angles which can 
guarantee diversity improvement and no error floors, and Hi) for a quantized receiver with perfect channel 
knowledge and minimum distance decoding, we analytically show that, among all admissible rotation 
angles, a good choice is one in which the transmitted vectors have equidistant projections along both the 
transmitted components. We then show that the square M^-QAM constellation rotated by 6^ = tan^^(l/Af) 
has equidistant projections. 

Further, we relax the perfect channel knowledge assumption, and propose novel training sequences and 
channel estimation scheme, which achieve an error probability performance close to that achieved with 
perfect channel knowledge. Through Monte-Carlo simulations we show that even with coarse analog-to- 

'Error probability performance is said to floor, if and only if it converges to a non-zero positive constant as the signal-to-noise ratio tends 
to infinity. 
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digital conversion, and short training sequences, the error performance with the estimated channel is similar 
to that with perfect channel knowledge. The main interesting result is that, even when the channel estimate 
is not perfect, an error probability performance exactly same as that with perfect channel estimate is 
achievable under some sufficiency conditions on the channel estimate and the number of quantization bits. 
These conditions are analytically derived, and shown to be satisfied by the proposed training/estimation 
scheme for some scenarios. Another interesting result is that, with sufficient number of quantization bits, 
the error performance never floors irrespective of the quality of the channel estimate. 

IL System model and Quantized Receiver 

We consider SISO block fading channels with single transmit and single receive antenna. The channel 
gains are assumed to be quasi-static for the coherence interval of the channel, and change to an independent 
realization in the next coherence interval. We further assume that the signaling bandwidth is much smaller 
than the coherence bandwidth of the channel (frequency flat fading), and therefore the channel frequency 
response is assumed to have constant magnitude and linear phase within the signalling bandwidth. Let 
the radio frequency band used for transmission be (/c — W/2, + W/2), where fc^Wis the carrier 
frequency and W is the signaling bandwidth. The complex channel frequency response is then given by 

H{f) = \h\e-^'--^, |/-/,|<^, (1) 

and zero elsewhere, i.e., scaling by \h\, and a delay of r seconds. The transmitted signal is given by 

xit) = cos(27r/et) + X? sm{27i f,t))g{t - kT), (2) 

k 

where 1/T is the rate at which information symbols are transmitted, and Xk = + jx^ is the k- 
th transmitted information symbol. We assume pulse shaping signals which result in no inter-symbol 
interference (ISI) (e.g., g{t) = s'mc{Wt)). Prior to the transmission of K information symbols, there is a 
training phase in which a known preamble sequence of P symbols is transmitted to enable carrier frequency 
synchronization in the receiver (i.e., enabling the phased locked loop (PLL) in the receiver to lock to the 
transmitter's local oscillator) and also for tuning the receiver gain. In this paper we assume the preamble to 
be a constant amplitude carrier obtained by setting x(, = A and = in Q. The received signal during 
the training phase of duration P symbols, is given by y{t) = A\h\ cos{2it fct—2n /ct) J2k=o 9i{'t~''~)~kT). 
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Fig. 1. Receiver analog front end (APE) . 

Figure [H shows the signal path of the analog front end of a typical heterodyne receiver flS]. Let the 
combined gain of the Low Noise Amplifier (LNA), Mixer (MXR) and Low Pass Filter (LPF) be denoted 
by qafe- In the training phase, after the PLL has locked, the LPF output (I Path) is given by u\t) = 
AgAFE\h\ Ylk=o~^ 9{{t — 'T) — kT)+n\t), where n\t) is the white Gaussian noise in the receiver (I path). 
The LPF output is digitized using a Nyquist rate sample & hold type analog-to-digital converter (ADC), 
as shown in Fig. [TJ Let the input dynamic range of the ADC be —cj2 to Cq/2. We also refer to cj2 as 
the clip level, since any input greater than Cq/2 would be limited to Cq/2. For optimum performance, it 
is desirable that the range of the input signal to the ADC matches with the ADC dynamic range (ADC 
range matching). Due to fading, the input level at the ADC may vary, and therefore a variable gain 
amplifier (VGA) is generally used to ensure ADC range matching. The gain of the VGA is controlled by 
the automatic gain control (AGC) module (HI. During the training phase, the AGC detects the peak of 
the signal u^{t) using a conventional analog peak detector whose output is given by 



(3) 



Let X denote the peak absolute value of the transmitted symbols, x\ and x^, during normal 
information transmission phase. During information transmission phase, ADC range matching (i.e.. 



Y = 9vGA9AFE\h\^) requires the VGA gain to be 



9VGA 



Cq A 1 



(4) 



2 X Vagc—pk 

Since the ratio A/X and Cg/2 are known a priori, this computation is done in the AGC using simple 
analog circuits [0. In the rest of the paper, we assume that this computation is perfect. 



During the information transmission phase, the PLL tracking loop is turned off and the VGA gain setting 
is frozen to the value given by (|4]). Therefore, during this phase, the ADC input signal (I path) is given 
by z'{t) = gvGA{gAFE\h\ Ek=p~' dait -r-kT)+ n\t)) = f Y.k=p~^ ^dit -r-kT)+ gvGAn\t). 
The ADC input (Q-path) is similar. Subsequently, without loss of generality, we assume an ADC with a 
normalized clip level of Cq/2 = 1. Assuming perfect timing synchronization (i.e., receiver can perfectly 
estimate r), the k-th output of the sample & hold circuit, at time t = r + kT is given by 

/ / Q Q 

\h\X ' ^ \h\X' ^ ^ 

where wl = {r + kT) / qafe and w2 = rfi{T + kT) / qafe are i.i.d. Gaussian random variables 
with variance denoted by cr^/2. Let the average transmit power be denoted by Pt = ]E[|xfcp]. Then 
the instantaneous signal to noise ratio (SNR) at the output of the sample & hold circuit is given by 
linst = -Prl^P/c"^- Assuming a Rayleigh fading model with h ~ CA/'(0, 1) (Complex Gaussian with mean 
zero and variance 1), the average signal to noise ratio (SNR) is given by 7 = E/i[7j„st] = -Pr/o"^- The 
output of the sample & hold circuit is then quantized by a 6-bit uniform quantizer Q, as shown in Fig. [B 
The quantizer is modeled by the function Qbit),t E M, which is given by 



Qb{t) 



+1, m > (2'-' - 1) 

-1, ^{t) < -{2'-' - 1) (6) 



^^1^^, otherwise 



(7) 



where [xj denotes the largest integer not greater than x. For a n-dimensional complex vector z = 
(zi,Z2, — - ,Zn), let Q6(z) denote the n-dimensional component-wise quantized version of z. That is, 
z = (21, 22, ■ ■ ■ , Zn) = Q6(z) implies that 



Qb{z!) , ~zf = Qb{zf) t = l,2,...,n. (8) 



The k-th quantized received symbol, = rl + jr^ is therefore given by 

ri = Qbisi) , r'^ = Q,is^) (9) 
where si and are the real and imaginary components of the fc-th sample & hold output symbol. 
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Fig. 2. Achieving modulation diversity by coding across n different chiannel realizations. 

Modulation diversity coding is illustrated in Fig. |2l Coding is performed across n > 1 information 
symbols resulting in n coded symbols/codeword. These n coded symbols are interleaved and then 
transmitted over n independent channel coherence intervals (realizations). At the receiver, the channel 
outputs during the n coherence intervals are buffered, followed by de-interleaving and detection. Suitable 
coding across n independent channel realizations results in an n-fold increase in the diversity of reception. 
In fading channels, codes designed using algebraic lattices can achieve modulation diversity, and are 
therefore employed to improve the diversity of reception [|6]|. With an unquantized receiver, it is known 
that lattice codes based on algebraic rotations can achieve full modulation diversity [TTll lfTOl . However, 
with quantized receivers, this is no longer true. In this paper we consider the case of n = 2. Let the 
information symbol vector be denoted by u = (^1,^2)^, where the information symbols ui and U2 are 
restricted to square M^-QAM signal set, though a generalization to non-square QAM is trivial. Let the set 
Sm = {-{M - 1), . . . , -1, 1, ■ ■ ■ , (M - 1)} denote the M-PAM signal set. Then, M^-QAM is denoted 
by the set S]^ = {w + jv \ w,v E Sm}- The information symbols are coded using a 2 x 2 rotation matrix 
G, resulting in the transmit vector x = {xi,X2)^ = Gu, where 



G 



(10) 



cos(6') sin(^) 
— sin(6') cos(6') 

Due to QAM symmetry, one can restrict the rotation angle in (fTOl) to [0,7r/4). The set of transmitted 



vectors X and the peak component value X are given by 

A' = < X I X = Gu, ui, U2 G Sif \ , X = max < max max(|x[|, |) S (11) 

•X.&X j=l,2 L 

Also, let the channel gain during the transmission of xi and X2 be denoted by \hi\ and |/i2|, respectively. 
We assume hi and /i2 to be i.i.d. CA/'(0, 1). Let r = (ri,r2)^ denote the quantized received vector, where 
ri = r{ + jr^ and r2 = + jr^ are the ADC outputs during the transmission of xi and X2, respectively. 
From dS]) and ^ it follows that 

With the above quantized receiver model, maximum likelihood decoding is no more given by the minimum 
distance decoder, and is rather complex. Nevertheless, due to its lower decoding complexity, we shall 
assume a minimum distance decoder taking r as its input, and the output (detected information symbols) 
given by 

Gun 2 



u = arg mm 



diag(|/ii|, \h2\)(r 



X 



Gu\ \L / /Gu 



argmin^^ |r-(^) ) D, ( r - ( ^ ) | , D, ^ diag(l, p^) (B) 



u&SljXSjj \ \ X J I \ \ X 

where p = |/i2|/|^i| is the channel gain ratio, and f, ||.|| denote Hermitian transpose and Euclidean norm 
respectively. In the subsequent sections [ni] and |IVl assuming perfect receiver knowledge of p, we show 
that even with the suboptimal minimum distance decoder in (fT3l) . we can avoid error floors and also 
achieve modulation diversity. Finally in section |Vj we relax the perfect channel knowledge assumption 
and present a practical training scheme to estimate p. We show that even with coarse quantization (i.e., 
small number of quantization bits - b), the proposed channel estimation scheme achieves error probability 
performance close to that achieved with perfect knowledge of p. 

III. Rotation Coding in Quantized Receiver 

In an unquantized receiver, at high SNR, the word error probability is minimized by choosing the 
transmit vectors such that the minimum product distance between any two vectors is maximized 0. 
There also exists algebraic rotations which guarantee a non-vanishing minimum product distance with 
increasing QAM size jH. In this paper, we study the error performance of these rotated constellations 
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with a quantized receiver and minimum distance decoding, and derive the conditions under which full 
modulation diversity can be achieved. 

In case of a quantized receiver, the sample & hold outputs (|5]), are quantized to the appropriate 
quantization box containing it. As an example. Fig. [3] illustrates the rotated 4-QAM constellation with 
6 = 20°. The dark filled squares represent the 4 possible values taken by the real component of the 
normalized transmit vector x^/X = {x{/X ^ x^X)'^. The projections of the 4 possible vectors onto 
the first component (horizontal) are marked with a cross. A b = 2-bit quantizer is used along both 
codeword components. The dashed horizontal and vertical lines represent the quantization boundaries 
along the 2 components. As an example, in Figl3]the real component of the sample & hold output vector 

= (s{, ^2)^ (marked with a star), is therefore quantized to = (r(, r^)^ (Note that there are totally 16 
different quantized outputs marked with empty circles). The quantization box corresponding to the output 

is shown in the figure as a square with solid lines. 

As the noise variance — )• 0, the sample & hold output s is almost the same as the normalized 
transmitted vector x/X. Therefore at sufficiently high SNR, if there exists two different transmit vectors 
X and y, such that Qb(x/X) and Qb(y/X) are identical, then it is obvious that the error probability 
performance would floor as SNR -> 00. This is because, at high SNR the quantizer output would be the 
same irrespective of whether x or y was transmitted, which makes it impossible for the the receiver to 
distinguish between the two transmit vectors leading to erroneous detection. More formally, two transmit 
vectors x and y are said to be distinguishable if and only if Qb(x/X) ^ Qb{y/X). Therefore, in order 
to avoid floors in the error probability performance, we propose the first code design criterion. 

Criterion I .■ A necessary and sufficient condition to avoid error floors with a quantized receiver, is 
that any two transmit vectors must be distinguishable. 

To achieve full modulation diversity, it is required that even under deep fading conditions in one 
component, any two transmit vectors x and y must still be distinguishable in the other component. 
This, therefore, implies that the projections of all the transmit vectors onto any one component must be 
distinguishable by the quantizer in that component. Therefore, we have the second criterion. 

Criterion II .• Given a b-bit quantized receiver, in order to achieve full modulation diversity, a necessary 
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Fig. 3. Signal space at the quantizer input with b — 2 (real component). Rotated 4-QAM (6 = 20°) depicted with dark filled squares. 

condition on the rotation angle 6 is that, any two distinct transmit vectors x and y satisfy 

Qb{x^/X)^Qbiy^/X), t = l,2. (14) 

With a rotated M^-QAM there are totally distinct projections onto any component, and therefore the 
minimum number of quantization bits required for the transmit vectors to be distinguishable along any 
component is at least [2 log2(M)] . Hence, in order to achieve full modulation diversity a straight forward 
lower bound on b ia^ 

b> r21og2(M)l. (15) 

Subsequently, we assume that for a given M, b is fixed to the lower bound value in (fTST l. We further note 
that, with a. b = [2 log2(M)]-bit quantizer. Criterion II is not satisfied by all rotation angle j^. 

With a b = [2 log2(M)]-bit quantizer, the set of angles (between and 7r/4) which result in 
distinguishable projections along both the codeword components will be referred to as the admissible 
angles (i.e., angles which satisfy Criterion II). For example, with 4- and 16-QAM, the admissible angles 
lie in the range (tan^^(l/5) 7r/4) and (11.3° 16.9°), respectively. With increasing M, the interval of 

^ \x] denotes the smallest integer not smaller than x. 

^For example, even though 6 = l/2tan^^(2) guarantees a rotation code having non-vanishing minimum product distance, with ab — 4-bit 
uniform quantizer and Af ^ = 16-QAM it does not satisfy Criterion II. 
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admissible angles reduces. With 256-QAM, the range of admissible angles is only (3.47° 3.68°). Another 
interesting fact is that, for M^-QAM, 9 = tan^^(l/A<f) is always in the set of admissible angles. Further, 
as M increases, tan~^(l/M) ± e are observed to be the only admissible angles. 

Apart from the fact that the chosen angle must have distinguishable projections, it can be analytically 
shown that for Af ^-QAM, any rotation angle for which the rotated constellation satisfies 

Qi(x/X) = x/X , X e A" (16) 

does indeed achieve a diversity order of 2 (i.e., full modulation diversity since n = 2), with a b = 
[2 log2(M)]-bit quantized receiver and minimum distance decoding given by (fT3l) (See Appendix lAl and 
Theorem lA.ll ). Subsequently, a rotated constellation which satisfies (fT6l) shall be referred as being matched 
to the quantizer. It is easy to see that a rotated M^-QAM constellation is matched to a 6 = 2[log2(M)]-bit 
uniform quantizer, if and only if, the projections of the transmit vectors are component-wise equidistant 
and distinguishable. 

Even with a mismatched rotated constellation having distinguishable projections (i.e., when the 
projections are not equidistant), full modulation diversity may be achieved, but then the error probability 
would be higher, since some transmit vectors would be closer to the edge of their quantization boxes 
(making it easier for noise to move the transmitted vector to another quantization box when received) 
(illustrated through Fig|TO]in Appendix E). Following along the same lines as the proof in Theorem lA.U 
it can be shown that mismatched constellations result in a higher error probability when compared to 
matched constellations. This therefore leads us to the third code construction criterion. 

Criterion III .• In order to minimize the error probability of a rotated M^-QAM constellation with 
a b = \2log2{M)'\-bit quantized receiver, the rotation angle must be such that the rotated M'^-QAM 
constellation is matched to the quantizer. 

IV. Rotated Constellation Design for Quantized Receiver 

In this section, we construct rotated M^-QAM constellations which satisfy Criterion III. We had earlier 
observed that, for A/^-QAM, a rotation by 6^ = tan^^(l/M) appeared to be always in the set of admissible 
angles. In fact, it can be shown analytically that a rotation by 6' = tan~^(l/M), satisfies Criterion III (See 
Theorem IB.1[ Appendix |B]). 
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For M^-QAM with 9 = tan~^(l/M), it can be shown that the minimum product distance of the code 
is 4M/(M2 + 1) 4/M for M > 1). On the other hand, a rotation angle of 6^ = l/2tan^i(2) is 
known to have a minimum product distance of at least A/^/E irrespective of the QAM size. Also, for 
any rotation angle the error performance with a quantized receiver is inferior to that with an unquantized 
receiver. Hence, with increasing M, the error performance of a quantized recerver with 9 = tan~^(l/M) is 
expected to be increasingly less power efficient than that of a unquantized receiver with 9 = 1/2 tan"^(2). 
With increasing M, the set of admissible angles appeared to be only tan~^(l/M) ±e and therefore, it can 
be argued that, the best possible error performance with a b = [2 log2(M)]-bit quantized receiver would 
have a loss in power efficiency when compared to an unquantized receiver. However, this appears to be 
the cost to achieve full modulation diversity in quantized receivers with limited precision. 

V. Imperfect Receiver knowledge of p 

In the previous sections, in order to achieve full modulation coding diversity, minimum distance decoding 
at the receiver assumed perfect knowledge of p. In this section, we relax this assumption and present novel 
techniques to estimate p accurately. It is expected that the error performance would degrade with imperfect 
receiver knowledge of p. Interestingly, in sub-section IV-AI we propose an optimality criterion, which if 
satisfied by the estimate of p, would guarantee no loss in the error probability performance of the minimum 
distance decoder with estimated p when compared to the error performance with perfect knowledge of p. 
Such an estimate would be referred to as an optimal estimate of p. We estimate p based on the quantized 
receiver outputs for a known transmitted sequence. We refer to this transmitted sequence as the p-training 
sequence. Any p-training sequence which results in an optimal estimate of p is subsequently referred to 
as an optimal p-training sequence. In sub-section IV-BI we present receiver control techniques required to 
estimate p. Maximum likelihood estimation of p based on the quantized receiver outputs of the p-training 
sequence is discussed in sub-section IV-CI 

Finally, in sub-section IV-D[ for M = 2 (rotated 4-QAM) we present an optimal p-training sequence 
which satisfies the optimality criterion introduced in sub-section IV-AI For M > 2, the length of p- 
training sequences which satisfy the optimality criterion is expected to be large resulting in too much 
training overhead and hence loss in effective throughput. Therefore, a novel design of short p-training 
sequences is proposed, which can achieve an error probability performance close to that achieved with 
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optimal p-training sequences. Such short p-training sequences have been referred to as 'good' training 
sequences. Also throughout this section, it is assumed that i) with rotated Af^-QAM, ab = \2 log2(M)]-bit 
uniform quantizer is employed, ii) a minimum distance decoder is used for detection, and ///) the rotated 
constellation satisfies Criterion III. 

A. Criterion for the optimal p estimate 

Since the rotation matrix G is real-valued, it is obvious that the minimum distance decoder in (fT3T ) 
separates into independent and identical minimum distance decoders for the real and imaginary components 
of the transmitted information symbol vectors, and therefore the error probability performance for both the 
real and imaginary components are also identical. Hence, we only analyze the optimality of the minimum 
distance decoder with imperfect p estimate, only for the / component. The minimum distance decoder 
with the estimated p, is also given by (fT3] ). but with p replaced by its estimate p. 

We are now interested in studying the conditions under which the error probability performance with p 
is exactly the same as the error probability performance assuming perfect receiver knowledge of p. Any 
estimate of p, which satisfies these conditions would be an optimal estimate in terms of achieving an 
error probability performance same as that achieved with perfect receiver knowledge of p. To simplify 
notations, for any received vector r, information symbol vectors u and v and any real C > 0, we define 




(17) 




(18) 



The detected information symbols can therefore be stated in terms of m(.) as 



arg min m{p^v^ ^vJ). 



(19) 



This then implies that, for any information symbol vector v 



DB(p,r^u^v^ 



)<0 



(20) 



With an estimated p, if for all information symbol vectors v G 5|, 



M 



DEipy^vi'y 



)<0 



(21) 
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then it is obvious that the output of the minimum distance decoder with estimated p is the same as the 
output of the minimum distance decoder with perfect knowledge of p. If (|2TI) holds for all information 
symbol vectors v G Slj, then along with (l20l) . it follows that 

De{p. v\ u', ^')De{^ r^ vO > (22) 

for all information symbol vectors v G Since could be any information symbol vector in iS|^ and 
De{C, r^, v^) = —De{C, v^, u^), it is easy to see that the output of the minimum distance decoder 
with estimated p would be the same as that with perfect knowledge of p if 

De{p, r', u', ^')De{p, r', , ^r') > (23) 

for all possible received vector r (finitely many due to receiver quantization) and all possible information 

symbol vectors u and v. We formally prove this observation in the following theorem. 
Theorem 5.1: For a given realization of p, and estimated p, if (|23l) is satisfied for all possible received 

vector r and all possible information symbol vectors u and v, then p is an optimal estimate of p. 

Proof. See Appendix O ■ 
We now analyze the condition set-forth in Theorem 15 . 1 1 regarding the optimal estimate of p. With each 

information symbol belonging to Af^-QAM, and h = [21og2(M)] we make the following definitions 

A r (ai - 02) I ^ ^ \ „2 A r 2 1 ^ 1 

Qm = { _ I Qi> Q2, Qs, ^4 e , Qs 7^ Q4} , Qm = |a I a G Qm , a > o| (24) 

where Sm^ is the M^-PAM signal set (see Section |Il| for PAM set definition). It is also noted that Sm^ 
is not the same as Sjj. As an example, with M = 2 and b = [21og2(M)] = 2, Sm'^ = {—3, —1, 1,3} 
and Vj^j = {—2, —4/3, —2/3, 0, 2/3, 4/3, 2}. Based on Theorem 15. 1[ the next theorem gives a useful 
sufficiency condition for the optimality of an estimate of p. 

Theorem 5.2: Consider a rotated Af^-QAM constellation matched to the quantizer. Let p be an estimate 
of p satisfying the following condition 

yieQlj : < / ^ < / 

V/ G Q+ : p' > / ^ p' > / (25) 

Then, p is an optimal estimate of p. 
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Proof : See Appendix |Dl ■ 
Let = {0, gi, g2, ■ ■ ■ , QLm}^ "^ith < gi < g2 < ■ ■ ■ < QLm- It is clear that the elements of the set 
Q\i partition the positive real line [0, oo) into (La/ + 1) intervals with [0, qi) and [gL^/, cso) being the first 
and the last interval respectively. The /c-th intermediate interval is given by [g^^i, qu), /c = 2, 3, ■ ■ ■ , Lm- 
For any finite set S = {si, S2, ■ " " ? with < si < S2 <■■■<■§„, let I{S) be the set of intervals 
induced by the set S. That is 

I{S) = |[si,S2), [S2,S3), ■■■ , [S„,CX))| (26) 

The sufficiency condition in (|25] ) can now be understood in terms of the (Lm + 1) intervals of the 
positive real line induced by the set The sufficiency condition basically states that, for an estimate 
of p to be optimal, it must belong to the same interval of X(Q|j) in which p lies. This therefore also 
implies that, for an estimate p to be optimal it is not necessary that p be exactly equal to p. 

B. Receiver control for estimating p 

In this section, we discuss receiver control techniques required for estimating p. In the proposed rotation 
coding scheme, coding is performed across n = 2 channel realizations. Using a 2-dimensional rotation 
matrix G, a pair of information symbols is transformed into a pair of coded output symbols. The first 
coded symbol in the pair is transmitted during channel realization 1 (with channel gain whereas 
the second coded symbol is transmitted during channel realization 2 (with channel gain |/i2|)- For both 
channel realizations, the preamble sequence used for tuning the VGA gain is the same as discussed in 
Section |Ill However, for channel realization 2, even before transmitting this preamble sequence, a known 
p-training sequence is transmitted for estimating p. During the transmission of the p-training sequence 
the analog gains Qafe and Qvga are set to the values programmed during the transmission of coded 
information symbols in channel realization 1, and therefore 

The p-training sequence is a sequence of / distinct positive valued symbols with each symbol being 
transmitted multiple times to average out the effect of receiver noise Q Subsequently, we shall denote an 

''At the receiver, depending upon thie repetition factor of the p-training symbols, the cut-off frequency of the LPF is appropriately reduced, 
which helps in noise reduction. 
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arbitrary p-training sequence by T. Let the k-th training symbol be given by c^, k = 1, 2, ■ ■ ■ , /. The / 
corresponding inputs to the sample and hold circuit are given by 

Sk = 9AFE9VGA\h2\Ck , k = 1,2,- ■ -l (28) 

Using (1271 ) in (|28l ), the sample and hold, and quantizer outputs during the transmission of the p-training 
sequence in channel realization 2 are given by 

rk = Qb{s,) =Qb{p^),k = 1,2,... I. (30) 

In channel realization 2, after all the p-training symbols are transmitted, the receiver estimates p based 
on the / observations {r^, k = 1, 2 ■ ■ ■ , /}. 

C. Estimation of p 

The / discrete outputs of the ADC ({ri,r2, ■ ■ -ri}) can be used to estimate p as follows. Given the / 
discrete outputs, the maximum likelihood estimate (MLE) of p is given by 

pML = argmaxP(ri,r2, • • -nip, {ck},l) (31) 

p>0 

where P(ri,r2, ■ ■ - nip, {ck},l) is the probability that the / outputs take the values {ri,r2, ■ • ■r/} for a 
given channel gain ratio p, and the p-training sequence {ck}. For the k-th training symbol c^, since 
rk = Qb{pck/ X), from ^ it must be true that 

'"fc — 2^ , fk + 2^ .p 

< P < , if rk < 1 (32) 



Ck/X Ck/X 



1 1 



2^-1 



< p < 00 , if Tfc = 1 



The inequality in (|32l) defines an interval of the positive real line, which we shall denote by £fc , k = 
1,2, . . . ,/. Therefore the / outputs would be {ri,r2, ■ ■ - n} if and only if p G £, where C = P\^k=i^k- 
Further, we would call C as the "ML interval" corresponding to the training sequence {c^} and the / 
outputs {ri, r2, ■ ■ ■ r/}. Also, given the / outputs, all the values of p in the interval C are equally probable. 
Let Csup = sup £, and Cinf = inf C, denote the supremum and infimum of the interval C. One possible 
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ML estimate of p, that we propose, is then given by 



P 




otherwise. 



sup 



(33) 



For a given p-training sequence T = {ci, C2, ■ ■ ■ , q}, and a set of corresponding outputs IZ = 
{ri, r2, ■ ■ ■ ri}, let C{T, TZ) C IR+ denote the ML interval. 

As an example, let us consider a 6=2-bit quantizer, and a training sequence {cfc} = 
{X/4, X/2, X, 2X, 4X}. Let the / = 5 corresponding output symbols of the quantizer be {ri = 1/3, r2 = 
1/3, ra = l,r4 = l,r^ = 1}. The intervals Ck corresponding to these 5 outputs are 



The ML interval i:({X/4, X/2, X, 2X, 4X}, {1/3, 1/3, 1, 1, 1}) = [2/3, 4/3) and hence p = 1. 

For a 6-bit quantizer and some fixed p-training sequence T = {ci,C2, ■ ■ ■ ,q} of / training symbols, 
it is clear that for each value of p G [0, oo), there is a corresponding output sequence 7^(p, T) = {r^ = 
Qb{p^) ? ^ = 1? 2, ■ ■ ■ , /}. We shall refer to each such possible output sequence as a feasible output 
sequence for the given training sequence. Note that even though the range of values of p is infinite, the 
number of distinct feasible output sequences is finite due to the finite length of the p-training sequence 
and the finite number (2^) of quantizer levels for a 6-bit uniform quantizer. Further, for each feasible 
output sequence TZ', there exists an ML interval C{T, TZ'). Also, it is trivially true that two ML intervals 
corresponding to two different feasible output sequences are disjoint (i.e., if TZ' ^ IZ" then C{T,TZ') fl 
C{T,TZ") = 0, where denotes the null set). In addition, since p G [0,oo), it follows that the ML 
intervals corresponding to all possible feasible output sequences, form a partition of the positive real line. 
This is summarized as follows. For any < pi, p2 < oo , pi 7^ P2 




- < p < 00 , £4 : - < 2p < 00 , £5 : - < 4p < 00. (34) 




£(r,7^(pl,r)) = £(r,7^(p2,r)) , 7^(pl,r) = 7^(p2,r) 



, 7^(pl,r)^7^(p2,r) 



(35) 



Also, 




(36) 
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Some more interesting properties are as follows. For any p-training sequence T, we have the following 
interesting properties. 

p e c{T,n{p,r)), < p < cx) 

Pi e C{T, Tl{p2, T)) =^ P2 ^ C{T, Tl{pi, T)) and vice- versa < pi, p2 < oo , pi 7^ P2 

p e c{r,n{p,r)) , o<p<oo (37) 

where p is the ML estimate of p given by (1331) . For any T, p' G jC{T, 7^(p, T)), if and only if TZ{p, T)) = 
TZ{p', T)). This equality is trivially satisfied for p' = p, which completes the proof for the first statement in 
(l37l) . For the second statement, we note that pi G C{T, Tl{p2, T)) if and only if Tl{pi, T)) = 7^(p2, T)). 
This then implies that p2 G C{T, Tl{pi, T)). The third statement follows from the fact that, for any p, the 
ML estimate p is always between the supremum and the infimum of the interval C{T,TZ{p,T)). 

Given a fixed p-training sequence T, from (l37l) it follows that both p and the proposed ML estimate p 
lie in the ML interval C{T, 7^(p, T)). In addition to this, if the ML interval corresponding to each feasible 
output sequence is a subset of some interval induced by the set Q^^, then for any / G Qj^j it follows 
that if / is greater than p then it is also greater than p, and similarly when / is smaller than p then it is 
also smaller than p. However this is precisely the sufficiency condition in Theorem 15.21 We can therefore 
conclude that with the proposed p estimation technique (see (|33])). a training sequence T results in an 
optimal estimate of p, if T satisfies the following conditions. 

VO<p<oo : /:(r,7^(p,r)) C /, for some/ gX(Q|J. (38) 

We next design optimal and near-optimal p-training sequences based on the criterion in (|38). 

D. Design of p-training sequence for estimating p 

We first show that with M = 2 and a 6 = 2-bit uniform quantizer it is possible to design a p-training 
sequence which satisfies (l38l) and is therefore optimal. 

Theorem 5.3: Consider a rotated constellation matched to the quantizer. Let M = 2, 6 = 2 and 
Oil = {0,gi,g2, ■ ■■qiM}^^ < qi < q2 < ■■■ < qiu- The following p-training sequence {ck} of length 
Lm with the proposed ML estimator (Section IV-CI l results in an optimal estimate of p. 
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Proof. See Appendix |El ■ 
For a general M > 2, it is challenging to design an optimal p-training sequence which satisfies the 
sufficiency condition in Theorem 15.21 Further, we conjecture that, just as with M = 2, for any M > 2 also, 
the length of optimal p-training sequences based on Theorem 15.21 would be proportional to Lm = \ Qm\- 
However, the cardinality of is a rapidly increasing function of M (e.g., = 29, and \ Q'^\ = 3939), 
which then implies that with increasing M a significant amount of communication bandwidth would be 
used up in the transmission of the p-training sequence, resulting in reduced overall throughput. Therefore, 
it is of practical interest to design p-training sequences which are short and which can still achieve an 
error performance comparable to that achieved with optimal p-training sequences. 

Towards designing such practical sequences, we observe that the average error performance with 
estimated p, would be sensitive to the amount of overlap between the ML intervals inducecjfl by the 
p-training sequence and the intervals induced by Q^. With short p-training sequences, the ML intervals 
induced by the p-training sequence would not coincide exactly with the intervals induced by 
Nevertheless, it may be possible to design short length p-training sequences for which some of the ML 
intervals belong to X(Q|^). With short p-training sequences, an interval induced by Qj^j, which is exactly 
the same as some ML interval induced by the p-training sequence, shall be referred to as "covered" by 
that p-training sequence. 

From Fig. |4l we try to gain more insights into the problem of designing shorter length p-training 
sequences for M = 2. We observe that the density of the intervals induced by Q'^j (depicted with cross 
'X' marks on the horizontal axis) is much more higher near the origin than farther away. Furthermore, 
from the p.d.f. of p = |/i2|/|^i | (Rayleigh faded hi and /i2), we observe that most of the probability mass 
is distributed near the origiro. Based on these observations and the sufficiency conditions in Theorem 
15. 2[ it can be argued that, to have an error performance comparable to that of an optimal p-training 
sequence, any short p-training sequence should aim to "cover" the intervals of X(Q|j) which are closer 
to the origin. This reasoning is supported by two facts. Firstly, with i.i.d. channels gains \hi\ and |/i2|, 
the probability of p taking large values is small, and hence large values of p are expected to have lesser 
contribution to the average error probability than smaller values of p. Secondly, when p >> 1, any error 

^These are basically the ML intervals corresponding to all feasible output sequences for the given p-training sequence. 
*In fact, for any other fading distribution also, it can be shown that P{p < 1) = P{p > 1) = 1/2. 
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Fig. 4. The intervals induced by the set Q'^j for AI — 2. ML intervals for the proposed p-training and estimation schemes. 

in the estimation of p is likely to have a lesser impact on the error performance compared to when 
p < 1. To see this, we note that for p >> 1, the ML estimate for any p-training sequence would be the 
infimum value of the ML interval corresponding to the all ones output sequence which would also be 
large i.e., p >> 1. Therefore for any two transmit vectors x = (xi,X2)^ = Gu and y = {yi,y2)^ = Gv, 
i:)£;(p,r^u^v^)L)i5(p,r^u^vO ^ pV^^i > 0, where rfa = (4 - xi/Xf - (r^ - yi/Xf. Using 
Theorem I5.1[ this then implies that, with high probability, the output of the minimum distance decoder 
with estimated p is the same as its output with perfect knowledge of p. 

Therefore, any short p-training sequence should aim to "cover" the intervals of I{Qlj) which are closer 
to the origin. With M = 2, in Fig. |4]a short p-training sequence of length / < Lm is designed in a way to 
"cover" only the intervals of X(Q^^) which are closer to origin. This is done by non-uniformly sampling^ 
out / distinct elements of the set Q^j, such that the intervals induced by these / elements coincide with 
most of the intervals of X(Q|j) which are near to origin. Let us denote this set of / elements as Qj^j^ny 
The corresponding short p-training sequence which has ML intervals coinciding exactly with the intervals 
induced by Q^j^.^-^ is then given by (|39| ). where the set Q\j is replaced by the set Q^j^.j^ and Lm is 
replaced by /. In Fig. HI for / = 9, the elements of one such Q^^i] depicted through 'triangles'. Short 

^we use the word "non-uniform" since the sampling is biased towards choosing more elements which are closer to the origin. 
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p-training sequences which achieve an error performance close to that achieved with optimal p-training 
sequences, would be subsequently referred to as 'good' p-training sequences. 

Even though 'non-uniform' sampling of is one possible method for designing 'good' p-training 
sequences, with increasing M, the number of ways in which 'non-uniform' sampling can be done, would 
also increase rapidly, thereby increasing the complexity of finding 'good' p-training sequences. Therefore 
for large M, a simpler strategy is required to search for 'good' p-training sequences. We next present 
a very simple and parameterizable short p-training sequence design, which results in 'good' p-training 
sequences. The A;-th symbol of the proposed p-training sequence is given by 



Cfe 



Xd^^-^) , A;= 1,2,--- ,/ (40) 



where / is the length of the training sequence, and d > 1 is the ratio between the consecutive p-training 
symbols. Let us denote this p-training sequence by Td- This design is based on the observation that for 
many 'non-uniformly' sampled subsets Q^/;) C it is possible to find a value of d, such that a p- 
training sequence designed using (|40] |. would have ML intervals "almost" same as the ML intervals of the 
p-training sequence designed using the 'non-uniformly' sampled subset. That is, for many non-uniformly 
sampled subsets Q^j^) C for any p > 0, there exists some d > \ and some e close to zero, such 
that 

|£(r,,7^(p,r,)) n£(r^7^(p,r))1 < e\c{T\n{pX))\ m 

where T' refers to the p-training sequence designed using the given non-uniformly sampled subset C 
Qm F*^^ M = 2, this fact is illustrated through Fig. SI where for the given 'non-uniformly' sampled 
subset of Q\j, i.e., Q(^^,) (shown with triangles), a p-training sequence designed using (|40l ) has ML 
intervals (shown with circles) almost coinciding with the intervals induced by Q^^^;). For a given M and 
/, the optimal d can be found at reasonable complexity, by minimizing the average error probability (as 
a function of d) using Monte-Carlo techniques. 
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Fig. 5. BER vs. SNR for a quantized receiver, b — 4, 16-QAM and perfect channel state information at receiver. 

VI. Simulation Results 

All error probabilities reported in this section have been averaged over the Rayleigh flat fading statistics 
of the channel. Also, the receiver is assumed to have perfect channel state information. In Fig. [51 we plot the 
average bit error rate/probability (BER), for rotated 16-QAM constellation (M=4) and a 6=4-bit quantized 
receiver. The following four important observations can be made in Fig. [51 i) with 9 = l/2tan"^(2) 
(which is known to achieve full modulation diversity in an unquantized receiver, but does not satisfy 
Criterion II), the BER performance with a quantized receiver fails to achieve full diversity (note the 
difference in slope at high SNR), which validates Criterion II, //) with 9 = tan~^(l/4), which results in 
equidistant projections, the quantized receiver achieves full modulation diversity with 6=4. Further, the 
quantized receiver performs only 1 dB away from an ideal unquantized receiver at a BER of 10^'', Hi) 
with a quantized receiver a rotation angle of 6* = 16° also appears to achieve full modulation diversity, 
but perform poor when compared to a matched rotated constellation with 9 = tan~^(l/4). This supports 
Criterion III, and iv) In Fig. [Hit is also observed that with 16-QAM rotated constellation (9 = tan~'^(l/4)), 
the error performance floors with 6=3 < 4 quantization bits, which validates code design Criterion I. 

It was discussed in Section [IVl that with increasing QAM size, a quantized receiver would be 

*For any real interval I, \I\ = (sup / — infJ) refers to the length of the interval and I'^ refers to the complementary set R — / (i.e., all 
real numbers which do not belong to /). 
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Fig. 6. BER comparison between quantized and unquantized receivers. 4-,16-,64-QAM. Perfect ciiannel state information at receiver. 

increasingly less power efficient when compared to an unquantized receiver. This fact is illustrated in 
Fig. [6l where the BER performance of both unquantized receiver with 6'=l/2 tan-^(2) and quantized 
receiver with 6'=tan^^(l/M) are plotted for A'P=4-,16- and 64-QAM and b=\2 log2(M)] . Perfect channel 
state information is assumed at the receiver. It is observed that for a fixed BER of 2 x 10^'', with 
increasing QAM size, the increase in signal power required by a quantized receiver is more than that for 
an unquantized receiver. An unquantized receiver requires 6.3 dB more transmit power when the QAM 
size is increased from 16 to 64. For the same increase in QAM size, a quantized receiver would require 
7.8 dB more transmit power. However, when the QAM size is increased from 4 to 16, the extra transmit 
power required for a fixed target BER of 2 x lO^'^ is roughly the same (about 7.7 dB) for both quantized 
and unquantized receiver. In FigJTTI (Appendix lA]). we report the BER performance of a rotated 16-QAM 
constellation for varying 9 and fixed SNR. It is observed that, the rotation angle 9 = tan~^(l/M) (which 
results in a matched constellation) achieves the minimum BER. This then supports code design Criterion 
III. 

In Fig. 171 the BER performance with minimum distance decoding and imperfect knowledge of p, is 
plotted as a function of 7 for M = 2. The rotation angle is 6^ = tan^^(l/2) and b = 2. We firstly make a 
note that, for a matched rotated constellation, the error probability performance with p being any arbitrary 
positive valued estimate of p, does not have error floors. This is because, in the absence of noise (i.e.. 
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Fig. 7. BER perfoimance with a quantized receiver {b = 2) and imperfect receiver knowledge of p. Rotated 4-QAM (M = 2). 

7 = oo), when a certain information symbol vector v is transmitted, with r as the quantized output vector, 
the detection metric of some information symbol vector u (i.e., m(p, r^, u^) and m(p, r*^, u*^) ) is equal 
to zero only for u = v, and is positive for all other possible information symbol vectors. Therefore, the 
detected information symbol vector u is the same as the transmitted vector v, resulting in zero probability 
of error. This argument is supported by the fact that in Fig. El a fixed estimate of p = 1, has no floors in 
its BER performance. 

In Fig. Ul we also observe that the BER performance of the optimal p-training sequence de- 
signed using Theorem 15.31 is same as the BER achieved with perfect knowledge of p. Fur- 
ther, a short (/ = 9) 'non-uniform' sampling based p-training sequence designed with Q^;-, = 
{1/9,1/5,1/4,4/9,5/8,1,5/3,8/3,4} C achieves a BER is close to that achieved by the optimal 
p-training sequence designed using Theorem 15. 3 1 (compare curves 3 and 4). Also, the BER performance of 
a 'uniform' sampling based p-training sequence design with Q^^^^-, = {1/9, 8/9, 8/5, 9/4, 3, 4, 5, 8, 9} (the 
induced intervals are almost uniformly distributed) is inferior to the BER performance achieved by a 'non- 
uniform' sampling based design (compare curves 4 and 5). Finally, it is observed that, the BER achieved 
with the p-training sequence designed using (|40|) (with = 1.57, / = 9) is similar to the BER achieved with 
perfect knowledge of p (compare curves 1 and 6). For higher order rotated M^-QAM, we proposed 'good' 
p-training sequences which are short and have near-optimal performance. We support this fact through 
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Fig. 8. BER performance with a quantized receiver (b = 4) and imperfect receiver knowledge of p. Rotated 16-QAM (M — 4). 

FiglHl where we plot the BER performance for a rotated 16-QAM constellation (9 = tan~^(l/4)), with a 
quantized receiver (b = 4) and minimum distance decoding at the receiver with imperfect knowledge of 
p. An estimate of p is computed based on the proposed ML estimation scheme discussed in Section 
The p-training sequence used for estimation is the same p-training sequence used in simulation curve 6 of 
FigjVl From curve 4 in the FigJS] it is observed that with a short p-training sequence of only 9 symbols, 
it is possible to achieve a BER performance comparable to the BER performance achieved with perfect 
knowledge of p (curve 2). This is also interesting since the same training sequence was also observed to 
be near-optimal with M = 2, 6 = 2. It therefore appears that the length of near-optimal/'good' p-training 
sequences does not increase significantly with increasing QAM size. One possible reason for this could 
be that with increasing QAM size, the quantizer resolution b also increases, which makes the estimate of 
p more reliable. 
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VII. Conclusions 

In this paper, we addressed the problem of achieving modulation diversity in fading channels with 
quantized receiver. For 2-dimensional modulation coding, through analysis we showed that in quantized 
receivers with perfect channel knowledge, algebraic rotations with equidistant projections can achieve 
modulation diversity with low complexity minimum distance decoding. We then relaxed the perfect channel 
knowledge assumption, and proposed novel channel training/estimation, which were shown to achieve an 
error probability performance similar to that achieved with perfect channel knowledge. 
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Appendix A 

A MATCHED ROTATED CONSTELLATION ACHIEVES FULL MODULATION DIVERSITY. 

The main result that we prove in this section is that, a matched rotated constellation is guaranteed to 
achieve full modulation diversity. For the sake of clarity, we formally define error probability and diversity. 
The average error probability (averaged over both the fading statistics and the transmitted vector) is given 
by 

where it is assumed that all transmit vectors are equiprobable and -Pe(x) denotes the average error 
probability (averaged over the fading realizations hi and /i2) when x is the transmitted vector and minimum 
distance decoding is performed on the quantized output (fT3T ). Further it is implicitly assumed that Pe(x) 
is also a function of the SNR 7. The diversity achieved is given by 

6 4 hm i^f^. (43) 

7-s>oo log 7 

Using the union bounding technique, P{'y) can be upper bounded as 

where Pe(x, y) is the average pairwise error probability of the event that the minimum distance decoder 
decodes in favor of y when x was actually transmitted. 
If it can be shown that the pairwise diversity 

. A logPe(x,y) 

d(x,y) = hm = 2 (45) 

7->oo log 7 

then from the union bound in (1441) it follows that the achievable diversity 5 is indeed equal to 2 (since 
with n = 2 and i.i.d. Rayleigh fading, the maximum achievable diversity is 2). The following theorem 
shows that the statement in (l45l) is indeed true for a matched rotated constellation. 

Theorem A.l: Consider a rotated M^-QAM constellation which is matched with the b = 2[log2(M)]- 
bit quantizer. A minimum distance decoder on the quantized output achieves 5(x, y) = 2 for any two 
distinct transmit vectors x and y. 

Proof. We need to show that the pairwise error probability between any two transmit vectors has a diversity 
order of 2. We only consider the real components, since the rotation matrix is real valued and therefore the 
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error probability of the imaginary component would be the same. This implies that the pairwise diversity 
in (|45] ) is the same as that between the real component of the transmitted vectors. 

Consider any two distinct transmit vectors whose real components are denoted by x = (xi,X2)^ and 
y = {yi^y2Y ■ since the projections along both the codeword components are equidistant, Q5(x/X) = 
x/X and the same is true for y. Further, it is always possible to find two other vectors x' = {x'l, Xg)^, and 
y' = {y'li 2/2)^ ^"'^h that, x'^ = xi, x'2 = y2, y'l = yi and y'2 = X2. It is also noted that, Qb{x'/X) = x'/X, 
and the same is true for y'. Effectively, x/X, y/X, x.' /X and y' /X form the four vertices of a rectangle 
(Fig. |9l). We next focus only on this rectangle, and consider the constellation with only 2 possible transmit 
vectors, namely x and y. Therefore, we now have a reduced system in which only 1 bit is communicated 
(since we have only 2 possible transmit vectors). Also, for the reduced system we design a non-uniform 
quantizer such that the four vertices of the rectangle correspond to the four possible quantized outputs 
(Fig. |9l). As shown in Fig. [9l the received signal space is partitioned into four quantization regions by a 
quantization boundary in each component (depicted in the figure by dashed lines, there is one quantization 
boundary along the horizontal component and another along the vertical component). We are specially 
interested in investigating the worst possible pairwise error probability when x is transmitted and the 
minimum distance decoder decodes in favor of y. Since we are interested in the worst case, we choose 
the two quantization boundaries (one for each component) of the reduced system to be the same as those 
two quantization boundaries in the original system which separate x and y along both the components in 
the original system and are also closest to x in the original system. 

If the error probability of the reduced system has a diversity order of 2, then since the original system 
has more quantized output points than the reduced system (in addition to the four points of the reduced 
system), it is obvious that the error probability between x and y in the original system would always 
be smaller than that in the reduced system. Therefore with this assumption, it suffices to show that the 
diversity order of Pe(x, y) in the reduced system is 2. For the reduced system there are two possible 
cases, either /) x/X and y/X are at the positions shown in Fig. |9l or ii) the positions of x/X and y/X 
are swapped with each other. Further for each case, the positions of x/X and y/X could also be swapped 
with the positions of x'/X and y'/X respectively. This therefore results in four possible scenarios, out of 
which we only discuss the scenario shown in Figj9l since the analysis and pairwise diversity for the other 
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Fig. 9. Reduced system model. 

three scenarios is same. Out of the four possible scenarios, let us denote the scenario shown in Fig. |9]by 
the anti-clockwise ordered tuple starting from the first quadrant i.e., (y',x, x',y). With this notation, the 
other three possible scenarios are given by the tuples (x', y, y', x), (y, x', x, y') and (x, y', y, x'). 

Lets consider that x is the transmitted vector in the reduced system. For the minimum distance decoder 
to decode in favor of y, the received symbol vector (r{, must lie in a quantization box whose quantized 
output is closer to y in terms of the Euclidean distance dE{., ■) defined below. The Euclidean distance 
between any two vectors u and v is defined as 



From Fig. |9]it is clear that when |/i2|(/U2 + -^2) > \hi\{^i + Ai), (i£;(y',x) < dE{y',y) and (i£;(x',y) < 
(i£;(x',x). That is, the quantized output y'/X is closer to x/X than to y/X, and x'/X is closer to y/X 
in terms of Euclidean distance. Therefore, the minimum distance decoder decodes in favor of y only 
when the received symbol vector lies in the quantization box corresponding to either x'/X or y/X. This 
in turn happens only when the noise along the second component is less than — |/i2|/U2- Therefore, when 
|^2|(/^2 + A2) > |/ii|(/ii + Ai), a decoding error occurs if and only if, the noise along the second component 
(^2) is less than — |/^2|Ai2• 
Similarly, when |/i2|(/U2 + A2) < \hi\{iii + Ai), a decoding error occurs if and only if, the noise along 




(46) 
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the first component (w{) is more than These two error events are summarized as 

El : ^2 < -|/i2|Ai2, and |/i2|(/i2 + A2) > + Ai) 

E2 : w( > , and |/i2|(/i2 + A2) < + Ai) 



(47) 



The error events can be written equivalently as 



El : \W2< — /i2niax |/i2| , \hi 



(ah + Ai) 
(/i2 + A2 



and \h2\in2 + A2) > \hi\{fii + Xi 



Eo 



J (/i2 + A2) / IL I IL I (/^l + -^1. 



(/ii + Ai) 
and |/i2|(/i2 + A2) < + Ai) 



(/i2 + A2 



(48) 



Let Fr{E) be used to denote the probability of some event E. Pe(x, y) can therefore be expressed asL 



Pr(Ei U E2 
Pr(Ei)+Pr(E2) 



(49) 



where the second statement follows from the fact that Ei and E2 have no common support (i.e., -Ei □£'2 
(p). We further define two more events 

F J ■ f (/^2 + A2)\ flA I lA + 

El : <^ - W2 > mm [ fi2, fJ'ij T—:]ma.x\\h2\,\ni\ 



(/^i + Ai) 



(Ai2 + A2) 



and |/i2|(/i2 + A2) > \hi\{fii + \i 



Eo 



r . / (/i2 + A2)\ hi, I I/, + 



and |/i2|(/i2 + A2) < \hi\{ni + Ai 



(50) 



Since Ei C -Ei and E2 C -E2 it follows that 



Pe(x,y) <E;,,,;,, Pr(Pi)+Pr(4 



(51) 



'with a slight abuse of notation, we still use Pe (x, y) to denote the pairwise error probability in the reduced system. 
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A careful inspection of the error events Ei and E2 reveals that in both cases, we are interested in an 
event when a zero mean and finite variance Gaussian random variable is greater than the same constant 
(i.e., min ^/i2, /ii j)- Therefore, using (l50t and (ISTT ). we have 



,2 



w 



o A . / (/i2 + A2) 



max 1/^2^ ! a I /ill 



A (/ii + Ai)^ 
(/i2 + A2)2' 



(52) 



Then 



Pr(w > /3v/^) = $ I 



2/322 



(53) 



where <l>(x) = ^= r°° e 2 rft. Using (l53l) and (l52l) . we have 



Pe(x, y) < 



$1 



2/32^ 



(54) 



With a simple algebraic manipulation on the right hand side term in (1541 ). it can be shown that 

1 



i^e(x,y)< 



F,, I -r^ I e 2 dx, 



(55) 



/27rJo ~V2/32 

where F^(x) = Pr(2 < x) is the cumulative density function of z. Under the Rayleigh fading assumption, 
\h\\' and |/i2p are i.i.d. exponential random variables with mean 1, and hence F^{x) = (1 — e^^)(l — e~Q ). 
Integrating the right hand side term in (|55T ). we obtain 



Pe(x,y)<-<|1 



1 



+ 



1 



(56) 



As 7 — )■ 00 the series expansion of (|56l ) yields 



P.(x.y) < 1^7- + 0(7-') 



< 



3P| 



^7"' + 0(7"') 



(57) 
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Too close to the edge of quantization box 




Fig. 10. Example of a rotated constellation that is mismatched with the quantizer (i.e., non-equidistant projections), b = 2, 9 — tt/12 and 
4-QAM. 



tioiQ. 



where o(.) denotes the little-o notatior 

Since the projections along both the codeword components are distinguishable and equidistant, /ii, ^2, 
Ai and A2 (and therefore a and (3) are non-zero positive quantities for any pair of transmit vectors (x,y). 
The asymptotic expressions in (|57T ) reveal that the pairwise error probability Pe(x, y) indeed achieves 
second order diversity (i.e., 5(x, y) = 2). ■ 
We next discuss the reason why a 'good' choice of the rotation angle must be one which results 
in a matched rotated constellation. Through FigHOl we illustrate the reasoning behind the fact that a 
rotated constellation matched with the quantizer (i.e., with equidistant projections) achieves a lower error 
probability when compared to a mismatched rotated constellation. We plot the four transmit vectors and 
the 16-possible quantized outputs, for a mismatched rotated 4-QAM constellation (9 = n/12), and a 
b = 2-bit quantizer. The projections of the transmit vectors onto the horizontal component (marked with 
a cross) are not equidistant (distances between the projections are 0.423, 1.154 and 0.423). Due to this, 
the transmit vectors are not at the centre of their respective quantization boxes. This then implies that a 
transmit vector is closer to one edge of its quantization box compared to the other edges. Therefore for 

'"Any function f{x) in a single variable x is said to be o{g{x)) i.e., f{x) = o{g{x)) if — > as 2: — > 0. 
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Fig. 11. BER of rotated 16-QAM at SNR 7 = 30 dB and varying 6 in the range of admissible angles (11.3°, 16.9°). b — 4 and minimum 
distance decoding. 

a given noise variance, when compared to a rotation code with equidistant projections (i.e., a matched 
rotated constellation), a transmit vector in a code with non-equidistant projections has a higher probability 
of being received in a different quantization box whose quantized output is closer to some other transmit 
vector. 

Through the simulation plots in Fig|TTl we experimentally support the validity of the fact that a rotated- 
QAM constellation matched with the quantizer has a lower error probability than a mismatched rotated 
constellation. The BER performance of a rotated 16-QAM modulation code is plotted for varying rotation 
angle 9. The SNR 7 is fixed to 30 dB and the rotation angle is constrained to lie within the range 
of admissible angles, i.e., (11.3°, 16.9°) for 16-QAM. Minimum distance decoding is performed on the 
quantized output (b = 4). It is observed that the best BER performance is attained at ^ = tan~^(l/4) which 
corresponds to a matched rotated constellation having component- wise equidistant projections. Hence, this 
observation clearly validates code design Criterion III. 

Appendix B 
Theorem IB. II 

Theorem B.l: With Af^-QAM information symbols, the rotation code with 6 = tan~^(l/M) satisfies 
Criterion III. 
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Proof. The corresponding rotation matrix G, with 6 = tan ^(1/M) is given by 

M 1 



G ' 



VM2+T 



-1 M 



(58) 



Since G is real-valued, it suffices to prove the equidistant projections property for the real component 
only. The information symbols u{ and U2 take values from the M-PAM signal set Sm- It is easy to 
show that with G given by (|58l) . there are distinct equidistant projections on both the codeword 
components. The pair (u{,U2) can take any of the values from the ordered sequence of values 

Si = {(-M + 1, -M + 1), (-M + 1, -M + 3), . . . (-M + 1, M - 1), (-M + 3, -M + 1), (-M + 
3,-M + 3),...(-Af + 3,M-1),...(M-1,-M+1),(M-1,-M + 3),...(M-1,M-1)}. It is easy 



to see that the value of the first component of the transmit vector, x{, increases in steps of 2/ V A'P + 1 
as (wfjMg) takes values sequentially from the set Slj. This then proves that the projections along the first 
component are indeed equidistant. 

Also, since the values in the Af-PAM signal set are symmetric around 0, it follows that the set of all 
values which the first component x{ takes, is same as the set of the values taken by the second 
component Hence the projections along the second component are also equidistant. ■ 

Appendix C 
Proof of Theorem 15.11 

For a given channel realization, let the received vector on the real component of the channel be r^. 
From (fT9l) . it is obvious that 

m(p, u^) < m(p, u^) (59) 



for any information symbol vector u. Using the definition of De{-, ■,■,■) in (|T8I) . (1591) can be written as 

DEip,r',u',u')<0 (60) 

for any information symbol vector u. 

Let us assume that, the estimate of p, p, satisfies (l23l) for all received vector r, and all information 
symbol vectors u and v. Hence, p must satisfy (l23l) for u = u and v = u. This, then implies that 



De{p, r^ u^ u')De{p, r^ u^ iiO > (61) 
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for any information symbol vector u. 
Combining (|60l ) and (|6T| ) we have 

DE{p,r',u',u') <0 (62) 
for any information symbol vector u. This however implies that 

m(p, r'^, u'^) < m(p, r'^, u^) (63) 

for any information symbol vector u. This then means that the output of the minimum distance decoder 
with the imperfect estimate p, is indeed the same as the output of the minimum distance decoder assuming 
perfect receiver knowledge of p. Since, this is true for any received vector r, it is obvious that the error 
probability performance of the minimum distance decoder with imperfect p estimate is the same as the 
error probability performance with perfect knowledge of p. Hence, any estimate of p, which satisfies (l23l) 
for all r,u and v, is indeed an optimal estimate. ■ 

Appendix D 
Proof of Theorem 15.21 

We basically show that if (l25l) is satisfied for all / G Q\j, then this implies that (l23l) is satisfied for 
all possible received vector r and information symbol vectors u and v. The optimality of p then follows 
from the application of Theorem 15.11 

For any given received vector r, and information symbol vectors u and v, assuming (|25] ), we would 
eventually show that the given r, u and v satisfy the sufficiency condition in (l23l) . For a given r, u and 
V, let di and d2 be defined as 

d. = (ri - p - (ri - yjr 

d. ^ {ri - - {ri - (64) 

where 

X = (xi, Xif = Gu , y = (yi, = Gv. (65) 

It is noted here that, since the transmitted vectors have the equidistant projections property, (r( — 
x)' (^1^ - x)' (^2 - x)' (^2 - x) ^ This then implies that, when da ^ 0, di/da G Qm- Further, 
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De{C, r^, v^) can now be written in terms of di and ^2 as follows 

D£(C,r^u^vO=cii + C'c^2. (66) 

We now distinguish between two important events, depending on whether = or (ii = or both 
are non-zero. We first consider the case when at least one among di and d2 is zero. If c/2 = 0, then using 
(|66l ) it is trivially true that De{p, r^, u^, v^)De{p, r^, u^, v^) = df > 0, thereby satisfying the sufficiency 
condition in (l23l) . Similarly, when di = 0, it is again trivially true that De{p, r^, u^, \^)De{p-, r^, u^, v^) = 
d^p^p^ > 0, thereby satisfying the sufficiency condition in (|23] ). 

When (^2 7^ and di ^ 0, we consider two cases depending on the sign of di and ^2. If both di 
and c/2 have the same sign, then sign{di + p^d-z) = sign{di) = sign{di + p^dz), and therefore the 
sufficiency condition in (l23l) is satisfied. Hence, the only case which remains to be handled is when 
sign{di) = —sign^dz). 

When sign{di) = —sign{d2), we have two possible situations, depending upon the sign of {di + p^dz). 
When (di + p^^) > and di > 0, ^2 < 0, it follows that p^ < {\di\/\d2\). Since {\di\/\d2\) > and 
(M1I/M2I) G Qti, from (125]) it follows that p^ < {\di\/\d2\). Since, di > 0, ^2 < 0, this then implies that 
{di + p'^d2) > 0, and therefore De{p, , , v^)De{p, r^, u^^, v^) = {di + p^c?2)('^i + P^c^2) > 0, thereby 
satisfying the sufficiency condition in (1231) . 

If {di + p'^dz) > and di < 0, ^2 > 0, it follows that > {\di\/\d2\). From (l25]), it then follows that 
p^ > (|(ii|/|(i2|). Since di < 0,^2 > 0, (rfi + p'^d2) > 0, thereby satisfying the sufficiency condition in 
(El. 

In a similar way, if (di + p^c/2) < 0, then also it can be shown that the sufficiency condition in (l23l) is 
satisfied. ■ 



Appendix E 
Proof of Theorem 15.31 



For the p-training sequence given by (1391 ) and the proposed ML estimation scheme, we would like to 
show that the ML interval of each feasible output sequence is exactly same as some interval induced by 
the set 
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From (|39l ) it is obvious that, since {g^} is an increasing sequence w.r.t. increasing k, {ck} is also an 
increasing sequence. Since the uniform quantizer function Qb{ ) given by Q is also a monotonically non- 
decreasing function of its argument, it follows that i) the output sequence {rk} is also a non-decreasing 
sequence , and ii) since {c^} is a positive valued sequence, it follows that {r^} is also a positive valued 
sequence. For a b = 2-bit quantizer, the only possible positive values that can be taken by the quantizer 
output are {1/3, 1}. These properties therefore imply that, the set of all feasible output sequences must 
be a subset of the set consisting of the following L^/ + 1 distinct sequences. 

n = r2 = ■ ■ ■ = tlj^i = 1, 
n = ■ ■ ■ = Tk = ^, rfc+i = ■ ■ ■ = rij^f = 1, k = 1,2,--- Lm - 1 

n = ra = ■ ■ ■ = Vl^^ = ^. (67) 

We will however show that, in fact, the set of all feasible output sequences is actually same as the set 
of sequences given by (l67l) . For each sequence in (l67l) . we derive the ML interval corresponding to it 
being a possible output sequence. We see that each ML interval is non-empty and coincides with exactly 
one of the intervals induced by Q'lj. This would then imply that, firstly, each output sequence in (|67|) is a 
feasible output sequence (since its ML interval is non-empty), and secondly, the estimate of p is optimal 
(since the ML interval of each feasible output sequence coincides with one of the intervals induced by 
Oil, and therefore from (|38T l it follows that p is optimal). 

We now derive the ML interval corresponding to each output sequence in (l67l) . Using (|32|) . the ML 
interval corresponding to the output sequence {ri = r2 = ■ ■ ■ = r^^^ = 1} is [2X/3ci, oo). From (|39| ). 
ci = 2X/?)qLj^., and therefore, the ML interval for the all ones output sequence is [qi^j, oo), which 
coincides with the last interval of Ql^. Similarly, the ML interval corresponding to the output sequence 
Ti = r2 = ■ ■ ■ = tlj^j = I is [0, 2X/3cl^j^). From (l39l) . cl^ = 2X/3gi and therefore, the ML interval 
corresponding to the output sequence ri = r2 = ■ ■ ■ = r^^^ = | [0)?i)' which coincides with the first 
interval induced by Q\.^. 

For the A;-th output sequence ri = ■ ■ ■ = = |,Tfe+i = ■ • • = tl^^ = 1, using (|32l) . the ML interval is 
given by [2X/3ck+i,2X/3ck) = [qL^-k, qiM-k+i)^ which corresponds to one of the intervals induced by 
Q^j. This then proves the optimality of the p estimate derived using the proposed ML estimation scheme 
with the p-training sequence given by (|39|) . ■ 



